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Abstract 

Background: The intestinal tract liarbours a complex and diverse microbial population that is important for health, 
yet has been poorly described in many species. This study explored the fecal microbiota of semi-free-ranging Wood 
bison {Bison bison othoboscoe). 

Results: A total of 2081936 16S rRNA (V4) sequences from 40 bison were evaluated. CatchAII analysis of richness 
predicted a mean of 10685 species per sample (range 5428-24764, SD 4136). Diversity was high, with an average 
inverse Simpson's index of 31.78 (SD 15.3, range 8.55-86.7). Twenty-one different phyla were identified; however, 
only Firmicutes and Proteobacteria, Actinobacteria accounted for >1% of sequences. Two distinct population clusters 
(Group A, n = 19 and Group B, n = 21) were evident based on both community membership and population structure. 
Group A had a significantly lower relative abundance of Actinobacteria (6.4 vs 1 1 .8%, P = 0.002), Chloroflexi 
(0.002 vs 0.013%, P = 0.014), Gemmatimonadetes (0.007 vs 0.15%, P = 0.038) and Proteobacteria (18.7 vs 42.5%, 
P= <0.0001) and a greater relative abundance of Firmicutes (70.9 vs 39.3%, P < 0.0001) than Group B. Within 
Group B, Alphaproteobacteria was the most common class of Proteobacteria (28% of all sequences), while 
Caulobacteraceae (18.5%), Pseudomonadaceae (3.5%), Hyphomicrobiaceae (3.5%), Alcaligenaceae (3.1%) and 
Xanthomonadaceae (2.6%) were the most abundant families. 

The twenty (3.1%) most abundant genera accounted for 71% of sequences. No operational taxon units (OTUs) were 
found in all samples at a relative abundance of 1% or greater. One OTU {Clostridium cluster XI) was present at 1% or 
more in all Group A samples, with two other Clostridium cluster XI OTUs in 18/19 (95%) samples. No OTUs were found 
at that abundance in all Group B sample, but an unclassified Lachnospiraceae was present in 20/21 (95%) and 
Clostridium cluster XI and Brevundimonas were found in 19 (90%) samples. 

Conclusions: The fecal microbiota of Wood bison is rich and diverse. The presence of two distinct populations not 
associated with housing, age or gender suggest that enterotypes, distinctly different microbial population 
compositions that can achieve the same ultimate function, might be present in bison, as has been suggested in 
humans. 
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Background 

The American bison (Bison bison) is a ruminant herbi- 
vore that is found in the wild in regions of western 
Canada and the United States and also raised in captivity 
as a food animal. Wood bison {Bison bison athabascae) 
are a geographic variant of Plains Bison {Bison bison bison) 
that are morphologically larger, but not genetically distinct 
and have historically inhabited areas in northern Canada 
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and Alaska. While they are raised commercially, little is 
known about the biology and health of bison compared to 
other domestic food animals. Some information can be 
inferred from other ruminants, but there are potentially 
important physiological and microbiological differences 
[1,2] that indicate the need for specific study. 

The gastrointestinal microbial population, the micro- 
biota, is receiving increasing study in many species be- 
cause of its role in health and various diseases, as well as its 
complex interaction with the body. Despite the known im- 
portance of the intestinal microbiota in many other species, 
there has been limited study in bison. Culture-based 
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methods have evaluated the fecal and ruminal microbiota 
of bison [1-4], but it is now clear that such approaches 
have significant inherent limitations for characterization 
of the microbiota. Culture-based methods provide a very 
superficial understanding of the intestinal microbiota 
because of its complexity, richness (number of different 
organisms), the fastidious nature of many microorganisms, 
the potential for profound culture bias (over-representation 
of bacteria that grow well under typical culture conditions) 
and because a large percentage of the microbiota is 
currently unculturable. While not without their own 
potential biases, culture-independent methods, parti- 
cularly next generation sequencing, have now become 
the standard for assessment of complex polymicrobial 
environments. 

Understanding the microbiota of normal animals can be 
important for assessing the role of the microbiota in 
disease, as well as the impact of various management 
practices (e.g. diet change, antimicrobial administration) 
on the microbiota. Accordingly, the objective of this 
study was to evaluate the fecal microbiota of members 
of a semi-free-ranging extensively managed wood bison 
herd using a culture-independent next generation sequence- 
based approach. 

Methods 

Study population and sample collection 

Fecal samples were collected from 40 semi free-ranging 
wood Bison, located in a 60 square kilometre fenced area 
within Elk Island National Park near Edmonton, Alberta, 
Canada. Bison are fenced behind a 2 metre high fence, 
but are essentially managed as wildlife within a national 
park. Animals were a mix of ages and sexes, but were 
fed the same hay diet for approximately one month 
prior to sampling. Bison originated from three separate 
pens in which animals were segregated by age and sex. 
Pen A housed 7 male calves (approx. 7 months old), 
pen B housed 20 adult females (>3 years of age) and 
pen C housed 13 adult males (>3 years of age). None of 
the animals had any history of antimicrobial exposure 
within the past two years. Fecal samples were collected 
per rectum during routine handling in a custom de- 
signed facility with alleyways and squeeze chute, which 
occurs biennially to allow management of surplus animals. 
Fecal samples were frozen at -20°C prior to analysis. This 
study was approved by the University of Saskatchewan 
Research Ethics Office (protocols 20130032 and 20130037). 

DNA extraction and quality control 

DNA was extracted from fecal samples using a com- 
mercial kit according to the manufacturers instruc- 
tions^. DNA quantity and quality were accessed by 
spectrophotometry^. 



16S rRNA gene amplification and sequencing 

The V4 region of the 16S rRNA gene was amplified 
using the primers forward S-D-Bact-0564-a-S-15 (5'- 
AYTGGGYDTAAAGNG-3') and reverse S-D-Bact- 
0785-b-A-18 (5'-TACNVGGGTATCTAATCC-3') [5] as 
has been previously described [6]. PCR products were 
evaluated by electrophoresis in 2% agarose gel and puri- 
fied with the Agencourt AMPure XP system.*^ Sequen- 
cing of the library pool was performed at the University 
of Guelphs Advanced Analysis Centre using an Illumina 
MiSeq^ and 2x250 chemistry. 

Bioinformatics 

The mothur package of algorithms (vl.32.1) was used 
for analysis [7]. Paired end reads were aligned. Se- 
quences were aligned with the Silva 16S rRNA refer- 
ence database (www.arb-silva.de) [8]. Sequences that 
were >244 bp or <239 bp in length, contained any 
ambiguous base calls or long runs (>8 bp) of holopoly- 
mers or did not align with the correct region were 
removed. Chimeras were identified using uchime [9] 
and eUminated. CatchAll was used to assess species 
richness [10]. Taxonomy was assigned using the RDP 
taxonomy database (http://rdp.cme.msu.edu/index.jsp). Se- 
quences were binned into operational taxon units (OTUs) 
at a 3% dissimilarity level. 

Subsampling was performed to normalize sequence 
numbers for further comparison. This consisted of ran- 
dom selection of a number of sequences from each sam- 
ple that corresponded to the lowest sequence abundance 
of all samples. Coverage was assessed using Goods 
coverage and visualized using rarefaction curves. Popula- 
tion diversity was described using the inverse Simpson s 
index and evenness was assessed using Shannons even- 
ness. Community membership was compared using the 
traditional Jaccard index, while community structure 
was assessed using the Yue & Clayton measure of dis- 
similarity. The core microbiota was assessed through 
identification of OTUs present in all samples at a mini- 
mum relative abundance of 1%. 

Visual assessment of dendrograms identified two dis- 
tinct groups, which were studied further. Relative abun- 
dances and alpha diversity indices were compared using 
t-test. The parsimony test was applied to the Jaccard and 
Yue & Clayton trees. Dissimilarity was also visualized 
using principal coordinate analysis (PCoA), and the statis- 
tical significance of the observed spatial separation was 
assessed using Analysis of Molecular Variance (AMOVA) 
and analysis of similarity (ANOSIM). A P value of <0.05 
was considered significant for all comparisons. 

Results 

A total of 2081936 V4 16S RNA gene sequences from 40 
bison passed all quality control filters. The number of 
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sequences per sample ranged from 17667 to 106284 
(mean 52049, SD 16896). CatchAU analysis of richness 
predicted a mean of 10685 species per sample (range 
5428-24764, SD 4136). Goods coverage ranged from 
0.943-0.976 (mean 0.962, SD 0.0097). Rarefaction curves 
based on subsampling of 17667 sequences per sample 
are displayed in Additional file 1: Figure SI. Diversity 
was high, with an average inverse Simpsons index of 
31.78 (SD 15.3, range 8.55-86.7). Shannons evenness values 
were an average of 0.60 (SD 0.046, range 0.477-0.683). 

Two distinct clusters (Group A, n = 19 and Group B, 
n = 21) were evident visually with both the Yue & Clayton 
(Figure 1) and Jaccard (data not shown) dendrograms. 
This was also evident visually using PCoA (Figure 2). 
Community membership and structure were signifi- 
cantly different based on parsimony test applied to both 
the Yue & Clayton [P = 0.034) and Jaccard [P = <0.001) 
trees. Significant differences in community structure 
were also observed with AMOVA {P < 0.001) and ana- 
lysis of similarity (ANOSIM) (P< 0.001). There was no 
difference in evenness between groups (0.644 vs 0.640, 
P = 0.78) and while there was not a significant difference 
in diversity (inverse Simpsons, 36.7 vs 27.3) it approached 
significance {P = 0.06). There was no association between 
pen-of-origin and Group {P = 0.53, Table 1). 



Twenty-one different phyla were identified (Table 1, 
Figure 3); however, only three (Firmicutes, Proteobac- 
teria, Actinobacteria) accounted for >1% of sequences 
overall. There were numerous significant differences in 
relative abundances of different phyla between groups 
(Table 2), with the most pronounced differences being 
with Firmicutes and Proteobacteria. Group A had a 
significantly lower relative abundance of Actinobacteria 
(6.4 vs 11.8%, P = 0.002), Chloroflexi (0.002 vs 0.013%, 
P = 0.014), Gemmatimonadetes (0.007 vs 0.15%, P = 
0.038) and Proteobacteria (18.7 vs 42.5%, P=<0.0001) 
and a greater relative abundance of Firmicutes (70.9 vs 
39.3%, P = <0.0001) than Group B. Six hundred fifty two 
genera were identified, with Clostridium cluster XI the 
most abundant overall (13%), followed by Brevundimonas 
(9.4%). The predominant genera in each group are 
presented in Table 3. 

Within Group B, Alphaproteobacteria was the most 
common class of Proteobacteria (28% of all sequences) 
followed by Gammaproteobacteria (8.0%) and Betapro- 
teobacteria (6.0%). Caulobacterales (18.5%), Rhizobales 
(7.4%), Burkholderiales (5.9%), Pseudomonadales (3.9%) 
and Xanthomonadales (1.6%) were the main orders, 
while Caulobacteraceae (18.5%), Pseudomonadaceae (3.5%), 
Hyphomicrobiaceae (3.5%), Alcaligenaceae (3.1%) and 
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Xanthomonadaceae (2.6%) were the most abundant 
families. 

A total of 652 different genera were identified. Most 
were rare, with the twenty (3.1%) most abundant genera 
accounting for 71% of total sequences. No OTUs were 
found in all samples at a relative abundance of 1% or 
greater. One {Clostridium cluster XI) was found at a rela- 
tive abundance of at least 1% in 38/40 (95%) samples, two 
{Clostridium cluster XI and an unclassified Lachnospira- 
ceae) in 37/40 (93%) samples and 4 (two Clostridium 
cluster XI, unclassified Lachnospiraceae, Brevundimonas) 
in 33/40 (83%) samples. 



Table 1 Pen classification of two different groups of fecal 
microbiota in bison 



Pen 


Description 


Group A 


Group B 


Total 


A 


Male calves 


2 


5 


7 


B 


Adult females 


10 


10 


20 


C 


Adult males 


7 


6 


13 


Total 




19 


21 


40 



When Group A was analysed separately, one OTU 
{Clostridium cluster XI) was present at 1% or more in all 
samples, with two other Clostridium cluster XI OTUs 
were present in 18/19 (95%) samples. No OTUs were 
found at that abundance in all Group B sample, but an 
unclassified Lachnospiraceae was present in 20/21 
(95%), while a Clostridium cluster XI and Brevundimo- 
nas were found in 19 (90%) samples. 

Discussion 

The fecal microbiota of bison contains a rich and diverse 
microbial population. It was rather surprising to find two 
distinct subpopulations since animals harbouring these 
subpopulations had been co-housed and fed the same hay- 
based diet for approximately one month. Typical factors as- 
sociated with gut microbiota alteration such as differences 
in diet, management, diet change, age, antimicrobial admin- 
istration and gastrointestinal disease were not present, there 
was no association with age or gender, and no other 
differences between these groups were evident. Reasons 
for this are unclear but one potential relates to the 
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Figure 3 Relative abundance of tiie predominant bacterial phyla in the feces of 40 bison. 



concept of enterotypes'. It has been proposed that the 
intestinal microbiota of most humans can be categorized 
into three enterotypes based on predominance of different 
bacterial groups, Bacteroides, Prevotella or Ruminococcus 
[11]. The basis behind this concept is that different popu- 
lation structures may be able to achieve similar function 
in the intestinal tract; however, it is unclear whether truly 
distinct enterotypes exist [12] and, if so, the relevance. 
Enterotypes have not been reported in other species, 
yet the presence of markedly different fecal microbio- 
tas amongst similar animals in this study suggests this 
phenomenon might occur in bison. This could be further 
explored by study of the functional genetic composition of 
the microbiota rather than relying solely on phylogeny. 

The clinical relevance of the group differences noted 
here is unclear. Consistent with the concept of entero- 
types, it is possible that there is adequate functional 
redundancy so that even profound differences in microbial 
populations can achieve the same functional results in the 
intestinal tract. However, study of the potential health im- 
pacts of these different microbial populations is warranted. 
This is perhaps most relevant in Group B, considering the 
high prevalence of Proteobacteria and relatively low abun- 
dance of Firmicutes. Proteobacteria tend to be of much 



lower abundance in the feces of healthy ruminants and 
other species compared to the relative abundance noted 
here [13-16], and increases in Proteobacteria are often 
associated with 'dysbiosis' and intestinal inflammation 
[17-19]. Increases in the relative abundance of Proteo- 
bacteria have been reported with various types of intes- 
tinal tract disease in diverse species such as humans, 
horses and dogs [13,14,18,20]. This study cannot evalu- 
ate the health impacts of Proteobacteria or other groups 
in bison, but the striking differences and concern about 
this phylum in other species raise some questions that 
should be addressed further. In contrast, the high abun- 
dance of Firmicutes in Group A is consistent with previ- 
ous reports in ruminants [15,16,21] and hindgut 
fermenters [14,22]. This phylum consists of a broad 
range of Gram positive bacteria with diverse functions 
and clinical relevance. It includes a large number of 
species with important roles in digestion and energy 
conversion and is expected to be present in large numbers. 

The dominance of a small number of phyla is consistent 
with other fecal microbiota studies, where a large number 
of phyla may be present with only a few accounting for 
more than 1% of sequences. Typically, Firmicutes, Bacteroi- 
detes, Proteobacteria, Actinobacteria and Verrucomicrobia 
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Table 2 Relative abundance of bacterial phyla found In 
the fecal microbiota of forty bison, divided into two 
different groups 

Phylum Group A (n = 1 9) Group B (n = 21 ) P value 

Acidobacteria 0.00008 0.00004 0.43 

Actinobacteria 6.4 1.2 0.002 

Bacteroidetes 0.49 0.65 0.26 

Chlamydiae 0.018 0.036 0.24 

Chloroflexl 0.0029 0.013 0.014 

Deferibacteres 0 0.00005 0.033 

Deinococcus-Thermus 0.00004 0.0023 0.40 

Elusimicrobia 0.000013 0.000008 0.72 

Fibrobacteres 0.016 0.018 0.80 

Fusobacteria 0.0029 0.0022 0.67 

Gemmatimonadetes 0.007 0.15 0.038 

Lentispaerae 0.0028 0.00005 0.24 

Plantomycetes 0.15 0.42 0.065 

Proteobacteria 18.7 42.5 <0.0001 

Spirochaetes 0.035 0.022 0.35 

Synergistetes 0.0001 0.0018 0.24 

Tenericutes 0.0015 0.0030 0.15 

Verrucomicrobia 0.54 1.3 0.056 

FIrmlcutes 70.9 39.3 <0.001 

SRI 0 0.000055 0.11 

TM7 0.42 0.42 0.97 

Unclassified 2.3 3.3 0.06 
Phyla that were significantly different area indicated in bold. 

are most common [14,23-26], as was the case here. The 
relative abundance of Bacteroidetes was low compared to 
studies of cattle [15,16,21]; however, care must be taken 
when making such comparisons and comparison of bison 
and cattle samples using the same laboratory and analytical 
methods is required to evaluate potential differences 
between those species. 

Three percent of genera accounted for 71% of sequences, 
yet no OTUs were found in all samples at a minimum 
relative abundance of 1%. Significant differences between 
the two Groups for most of the main genera may be one 
explanation of the lack of an identified core microbiota be- 
tween all bison. One Clostridium cluster XI OUT was the 
closest to a true core component, being found in all Group 
A samples and 90% of Group B samples as 1% of more of 
the sequences. Clostridium cluster XI has been reported to 
be a common constituent of the fecal microbiota of a di- 
verse range of species, including dogs, cats, pigs, grizzly 
bears and humans [27-30], although it has received limited 
investigation. It has been suggested that Clostridium cluster 
XI levels may correlate with dietary protein in carnivores; 
[29], however, this group also includes many saccharolytic 



organisms [31], something that might account for the rela- 
tively high abundance in hay-fed bison. It is likely that this 
is a very broad group from a functional standpoint, with 
individual members that can play different roles in diverse 
gastrointestinal environments. 

While this study can define microbial evenness and 
diversity, ideal values are not known so interpretation of 
these results in a clinical context is difficult. Greater diver- 
sity likely provides an added degree of functional redun- 
dancy, whereby there is a functional reserve' capacity in 
the microbiota to adapt to external influences. Yet, opti- 
mal diversity is unknown and some degree of unevenness 
(greater representation of some members of the micro- 
biota) is expected, since some bacteria play more im- 
portant roles than others (e.g. ceUulolytic organisms 
should probably be present at greater abundances than 
many other species in a herbivore). 

Comparison with previous culture-dependent study 
highlights the different results that can be obtained. One 
study of 96 bison samples involved a comprehensive 
culture-based approach yet only identified 19 'major' 
genera and a small but unclear number of lesser genera 
[4]. In contrast, 228 genera were identified in the current 
study, with species estimates exceeding 10000 per sample. 
Relative abundance was not evaluated in that study, but 
none of the ten most prevalence genera in the culture- 
dependent study were amongst the most abundant genera 
identified here (Table 3). Rather than true differences in 
populations, these contrasts likely reflect the difficulty 
growing many of the common genera identified here by 
sequence-based methods and overgrowth-bias of poten- 
tially uncommon genera that are adept at growing under 
conventional culture conditions. 

Fecal analyses must always be evaluated with an un- 
derstanding that there are differences in the microbiota 
of different parts of the intestinal tract. This may be 
particularly important in bison, as the proximal location 
and markedly different environment of the rumen mean 
that feces do not closely reflect the state of the ruminal 
microbiota [16]. Yet, while these differences may limit 
the sensitivity of the use of feces for detection of alter- 
ations in proximal compartments, it is reasonable to 
assume that differences in the fecal microbiota indicate 
differences in at least some proximal locations, and that 
the intestinal (post-ruminal) microbiota is also important. 
Comparative study of the entire gastrointestinal tract can 
be informative [16] but requires either surgical collection 
of samples or collection immediately after euthanasia, 
both of which carry inherent logistical, cost and ethical 
issues. Concurrent study of the rumen and feces is more 
practical because rumen fluid collection is more feasible 
(although certainly not without challenges in bison). 
Future study comparing the rumen microbiota between 
different bison, and comparing the rumen and fecal 
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Table 3 Relative abundance of predominant genera in the fecal microbiota of 40 bison 



Genus 


Total 


Group A 


Group B 


P value 


Clostridium cluster XI 


13.1% 


1 7.5% 


8.3% 


0.004 


Brevundimonas 


9.5% 


4.3% 


15.3% 


<0.001 


Unclassified Clostridiaceae 


8.9% 


10.1% 


7.4% 


0.013 


Unclassified Lachnospiraceae 


6.6% 


7.8% 


5.1% 


0.023 


Sedimentibacter 


4.2% 


5.9% 


1 .9% 


0.001 


Pseudomonas 


3.4% 


3.4% 


2.5% 


0.56 


U nclassified Peptostreptococcaceae 


2.9% 


3.9% 


1 .7% 


0.014 


Unclassified bacterium 


2.8% 


2.4% 


3.3% 


0.055 


Tissierella 


2.4% 


3.0% 


1 .9% 


0.16 


Devosio 


2.2% 


1.1% 


3.3% 


0.0007 


Unclassified Clostridia 


1.7% 


2.0% 


1 .5% 


0.065 


Unclassified Comamonadaceae 


1.6% 


1 .3% 


1 .7% 


0.55 


Nocordioidss 


\ .jyo 


1 . 1 70 


Z.DVO 


U.Uzo 


Plienylobocterium 


1.5% 


0.5% 


2.9% 


0.013 


Pusillimonas 


1.5% 


1 .0% 


1 .9% 


0.006 


Unclassified Ruminococcaceae 


1.4% 


1 .9% 


0.8% 


0.002 


Sporomuso 


1.4% 


2.5% 


0.7% 


0.008 


Aeromicrobium 


1.3% 


0.9% 


1 .8% 


0.18 


Cornobocterium 


1.3% 


1 .5% 


0.9% 


0.076 


Unclassified Bacillales 


1 .2% 


1 .5% 


1 .0% 


0.020 



microbiota would be useful Further, this study also 
only involved one population of bison, so extrapolation 
to other populations of wild or captive bison, particu- 
larly those under different management and from 
different regions, must be done with care. Study of wild 
and captive bison populations from different regions 
would be useful to compare with the results obtained 
here and to help determine the 'normal' microbiota and 
its variations. 

Conclusion 

Two distinct microbial populations were identified in 
a group of semi-free-ranging bison, with no apparent 
explanation for the differences. Numerous bacterial 
phyla were identified, yet a small number accounted for 
the vast majority of the microbiota. Further study of the 
potential that distinct enterotypes may exist in the bison 
fecal microbiota, and if so, the clinical implications of 
such as phenomenon, is indicated to further explore the 
potential impact of this on bison health. 

End notes 

^E.Z.N.A. Stool DNA Kit, Omega Bio-Tek Inc., 
Doraville, GA, USA. 
^NanoDrop, Roche, Mississauga, Canada. 
'^Beckman Coulter Inc, Mississauga, Ontario, Canada. 
*^Illumina, San Diego, CA, USA. 
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